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Abstract 

We propose an analytical scattering theory in spectral domain to model the electromagnetic 
(EM) fields of a gyrotropic sphere in terms of the eigen-functions and their associated spectral 
eigenvalues/coefficients in a recursive integral form. Applying the continuous boundary conditions 
of electromagnetic fields on the surface between the free space and gyrotropic sphere, the spectral 
coefficients of transmitted fields inside the gyrotropic sphere and the scattered fields in the isotropic 
host medium can be obtained exactly by expanding spherical vector wave eigenfunctions. Numerical 
results are provided for some representative cases, which are compared to the results from adaptive 
integral method (AIM). Good agreement demonstrates the validity of the proposed analytical 
scattering theory for gyrotropic spheres in spectral domain using Fourier transform. 

PACS numbers: 
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I. INTRODUCTION 



Electromagnetic scattering of anisotropic media have attracted more and more attention 
for their wide apphcations in the past decades, such as radar cross section (RCS) compu- 
tation of perfect electric conductor (PEC) targets coated with complex material, radome 
design, and interaction of light/wave with biological media and metamaterialsi""— . 

Based on the plane wave expansion in terms of spherical vector wave functions in isotropic 
mediumi^, the scattering by a uniaxial sphere and a sphere of uniaxial left-handed materials 
have been derivedi^>ii. More recently, the scattering of a gyromagnetic sphere has been 
investigated in the expansion in spatial domain^^. Moreover, the theory is only working for 
the case having gyrotropic permeability and scalar permittivity. If both permittivity and 
permeability are gyrotropic matrices, the interplay between the extra three parameters in 
the gyrotropic permittivity will make that approach^ too tedious and insufficient to model 
the scattering properties. This motivates our work in spectral domain instead of spatial 
domain. A most general gyrotropic sphere is considered, and since the existing method 
has drawbacks in the analysis of scatterings, a novel approach has to be developed. The 
analytical method, which can be readily implemented by programming, has its academic 
and practical significance in contrast to purely numerical solutions from FDTD, FEM or 
others. 

In view of this, we propose a distinguished method based on Fourier transform, and thus 
the spectral-domain analysis of the scattering by a general gyrotropic sphere in terms of 
spherical functions wave functions is investigated. This method has distinguished features: 
(1) it can straightforwardly be employed to describe the light wave interaction with particles 
and objects with gyrotropic permittivity and permeability; (2) the material constitution 
is very complex and general (both e and Jl are gyrotropic tensors), so all those existing 
scattering theorems are just its sub-cases, e.g., uniaxial, plasma, anisotropic, gyromagnetic, 
etc.; (3) it directly solves for the eigen-problems in spectral domain by Fourier transform, 
which simplifies the formulation in spatial domaiit^. 

To obtain the solution of vector wave functions in gyrotropic anisotropic media, we start 
from the vector wave equation in a source-free gyrotropic anisotropic medium. Taking the 
Fourier transform of the electric field and substituting it into the vector wave equation of 
the electric field, we obtain the characteristic equation. Solving this equation, the eigen- 
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values and corresponding vector wave eigenfunctions can be yielded. Then, electromagnetic 
fields inside and outside the gyrotropic anisotropic sphere can be expressed based on the 
eigenvalues and eigenfunctions. Those unknown scattering coefficients can be analytically 
determined from applying the continuous boundary conditions on the surface of the gy- 
rotropic anisotropic sphere, where orthogonality relations of the Legendre polynomials are 
employed. Numerical results are obtained to gain more physical insight into this problem. 
After the results were validated by comparison with the existing data, some new results are 
computed and discussed. 

In the subsequent analysis, a time dependence of the form exp{—iujt) is assumed for the 
electromagnetic field quantities but is suppressed throughout the treatment. 



II. ANALYTICAL FORMULATION 



The permittivity and permeability tensors of the gyrotropic anisotropic sphere shown in 
Fig. [T] are characterized by the following two matrices 

ei -ie2 
eg 



Hi —ifi2 
fis 



(1) 



The parameters are defined in Cartesian coordinates. The £^-field vector wave equa- 
tion can be obtained by substituting the above constitutive relations into the source-free 
Maxwell's equations^i, i.e.. 



V X [p;-^ ■ V X E{r)] - io^e ■ E{r) = 0. 
The solution to ([2]) can be obtained by the following Fourier transform: 



(2) 



E{r) 



dkr. 



E{k)e'^-''dh 



(3) 



where the wave number is denoted by fc = k^x + kyy + kzZ, and the space vector is identified 
as r = XX + yy + zz, with x, y, z being the unit vectors in Cartesian coordinates. By 




Free space 



FIG. 1: Geometry for the EM scattering of a plane wave by an gyrotropic anisotropic sphere. 



substituting (3) into (2), the wave equation can be transformed into 



where 



K{k) 



/OO poo POO 

dh / dky / K{k) ■ E{k)e^'^-^ dk, = 
-OO J —OO J —OO 



-bikj - b^k^ + ai h2,kxky - ih2kl - ia2 hik^k^ + ih2kyk^ 
bskxky + ib2k'^ + ia2 -bik'^ - bak^ + ai bikyk^ - ib2kxkz 
bik^kz - ib2kykz b^kyk^ + ib2kxkz -hik^ + kl) + 03 



with 



(4) 



(5) 



bi 

b2 
h 



/ii 



At2 



(6) 



In order to get nontrivial solutions of E{k), the following characteristic equation has to be 
satisfied: 

Det \K{k)\ = 0. (7) 



It can be explicitly rewritten as 

A{ekAk)k^ - B{dkAk)k'' + C = 0, 



(8) 



where 



C 



[6163 sin^ 6k + (&i — 62)^05^^^^] X [c^i sin^ 9^ + 03 cos^ 6k\ 
\bi{a\ — a\) + 630103] sin^ Ok + 203(6101 + 62^2) cos^ 6^. 
a-i{al - al) 



(9) 



with 



Ok 



k^ + ky + k^, 



tan-\^kl + kl/h), 
4>k = taxi'^ (ky/k^j,). 



(10) 



Equation ([8]) is a biquadratic equation with the following four roots of k^ (where £ = 1, 2, 3, 
or 4) for the radial wave vectors: 



k'^ 



k^ 



B + - AAC 
2A ' 



B - VB^ - AAC 
2A ■ 



(11) 



So the corresponding £^-field eigenvectors can be obtained from Eq. ([5]) and are given as 
follows 



k, <Pk 



F^l{9k, 4>k)x + F^(9k, h)y + F^\{9k, h)z 



f,{9kAk), (12) 



where and subsequently, g = 1, 2, 3, or 4; and 



Tpe Al A2 , 

= --^sin<pk + — cos0fc, 



qx 



qy 

qz 



—cos(pk + —sm(pk, 



(13) 



with 



Al 
A2 
A 



z(6i02 + b2ai)kgSin9kCos9k 

[bib3k'^sin'^9k + (61 — 62)^9^05^^^^] k'^sin9kCos9k — (6iai + 620-2) 

ai){hik'lcos^9k + h^k'lsiv?9k — cli) 



-{b2kqCos 9k + 02) + {bikgCos 9k 
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(14) 



With those obtained eigenvalues and their associated formulas, the E-fie\d in Eq. ([3]) is then 
given as follows 

E{r) = V / / Fl{ek,(l)k)Mek,(i)k)e''^''-''klsm9kd9kd(j)k (15) 

where 

kq = kq sin 6k cos + fc^ sin sin (l)ky + fcg cos O^z, 

and fq{9k,(t>k) denotes the unknown angular spectrum amplitude. Equation f[T^ is also 
known as the eigen plane wave spectrum representation of the electric field in homogeneous 
gyrotropic anisotropic medium. From ([3]), it is evident that the integration over the radial 
wave- vector component is reduced to a summation of four terms corresponding to the roots 
of dH]), which are the only permissible solutions. The symmetric roots, i.e., k = —kg oik = kq 
{q = 1,2) are taken into account automatically as 6 spans from to vr while spans from 
to 2tt. Physically, we need to sum up for only two of the four components, namely, ki and 

It is noted that the unknown angular spectrum amplitude fq{Ok, (pk) is a periodic function 
with respect to 6^ and 0^. Therefore we can use surface harmonics of the first kind to expand 

the fq{9k,(pk) 

WkAk) = Yl G'„,„,,P„^'(cosefc)e^™'^'= (16) 

m' ,n' 

where P^{x) denotes the associated Legendre function, n' is summed from to +oo, and 
m' is summed from —n' to n'. Substituting (fT6!) to (fT5|) . we obtain 

EirhYl E ^-'-'9 r r Fliek,(l)k)P:^' icosek)e'"''^e'^^''klsmekdekd(l)k. (17) 

q=l m',n' ^ "^^ 

This specific form of fll7p suggests the use of the well-known identity^"^ 



e^^-^= V^"(2n + l)Ukr) V i^i^^p-(cos^,)Pr(cos^)e^^ 
n-m)!_, im(^-^,) 



E f^^^r(cos^,)Pr(cos^)e 

m=l ^ ' 



(18) 



After substituting ( !T8|) into (fT7|) . we obtain the solution of E{r) for homogeneous gyrotropic 
anisotropic media. In order to have a compact and explicit solution to the scattering of a 



gyrotropic anisotropic sphere, it is necessary to introduce the spherical vector wave functions 
as follows^^ 



mn 



Pn (cos 6)^ rfP^(cOS^)- 



N, 



(0 



nin 



+ 1) 



Zn\kr 
kr 



sin 9 dO 



inuf) 



-.(0 



1 d{rz^\kr)) I dP^ {cos 9)-^ ^ .^P^ {cos 9) 
kr dr d9 sin 9 

rrfPr(cos^)- . P„"^(cos^)- 
" ^ -6 + im " . ^ ' cf) 



irncj) 



d9 " ' sin^ ^ " f ^^^^ 

where zi^{x) (where / = 1, 2, 3, or 4) denotes an appropriate kind of spherical Bessel 
functions, that is, jVi, 

or hn'' , respectively. Because of the complete property of the 
vector wave functions given in Eq. f lT9|) . we have the following expression 



ikaV 



^mnq {^kj-^"-^ mn 



M^^Sr,k,) + B^^^^{9k)N\^^{r,k,) 



+C:nn,{Ok)L^l{r,k,) 



(20) 



where n is summed from to +oo while m is summed from —n to n, and k is pointing in 
the (6'fc, direction while r is pointing in the {9, 0) direction in the spherical coordinates. 



The other inter-parameters. A, 



mnq \ 



B:,^^{9,) and C: 



mnq\ 



are provided in Appendix A. 



Substituting (120|) into (fT7|) . and integrating with respect to 0^, we end up with 



g=l 771,71 7l' 



mn'q 



A^,,(^,)M(^i(r, k,) + P^,,(^fc)ArW A;,) 



P„";^(cos^fc)fc„^in^fcrf0fc. 



(21) 



Equation f l2T]) is the eigenfunction representation of the £J-field in gyrotropic anisotropic 
media. The J/-field eigenvectors can be derived from ^^-field eigenvectors shown in Eqs. (j8])- 
f fTTj) by using the source-free Maxwell's equations in the spectral domain. Because the 
equations of iJ-field are very similar to those of £^-field, we only give the relation between 
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1?- field eigenvectors (i.e., F^) and ^^-field eigenvectors (i.e., F'i) in Cartesian coordinates 



UJ 



hi ib2 
-ib2 hi 
63 







sin 9k sin 0^ 



cos Ok 

' sin 9k sin 0^ sin 9^ cos 0^ 



— cos 9k 

— sin cos (pk 







(22) 



where g = 1, 2. 

From the result shown in (pTj) . it can be seen that the solutions to the source- free 
Maxwell's equations for the gyrotropic anisotropic medium are expanded in terms of the 
first kind of spherical vector functions. Because all spherical Bessel functions of different 
kinds satisfy the same differential equation and the same recursive relations, we can use 
the field expressions given in Eq. (1211) to analyze scattering and radiation by the stacked 
structure of the gyrotropic anisotropic media. 

Assume that the electric field of an incident plane wave is given by = xEoe'''°\ The 
incident EM fields (designated by the superscript inc) can be expanded by an infinite series 
of spherical vector wave functions for an isotropic medium as followsi^: 



(23) 



where 



.n+i 2n + l 

2n{n + 1)' 
,n+i 2n+l 

' 2 ' 
2n + 1 

2n{n + 1)' 
2n + l 



m 



m 



m 



m 



I s = l 
sj^l 



(24) 



According to the radiation condition of an outgoing wave and asymptotic behavior of spher- 
ical Bessel functions, only hn^ should be retained in the radial function, therefore the scat- 
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tering fields (designated by the superscript s) are expanded as 



mn 



(25) 



where and (with n being from to +00 and m being from —n to n) are unknown 
coefficients, and k^ = a;(eo/io)^''^) and /io denote the wave number, permittivity and 
permeabihty in free space, respectively. 

The expressions of EM fields inside the gyrotropic anisotropic sphere are given in Eq. (|2T|) , 
and the continuity of the tangential EM field components at r = a yields 

2 00 „^ 

2TTGmn'q / QmnqP^{cOS 6k)kl siu OkdOk 



where 







q=l n'=0 

= Eq [5m,l + Sm~l] 0. 
2 oo 



X 

mn 



{koay 



q=l n'=0 



RmnqPn' {cos Ok)kl siu 6kd9k 







^0 [5m,l + 5m,-l] K 



(koa) 



Q 



mnq 



X d 



+a 



ko 

h 3n{kqr 



mnq 



h^nKhr) 



R. 



mnq 



\ ko "'"'^ kor dr 
1 d 



k rdr ^^^'^^^'^^^^ 



~ mnq 



q' 

jn{kqr) 



h^n\kor) 



(26) 



(27) 
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The scattering coefficients, i.e., and are thus expressed as 

■ oo 2 



h'n\kQa) 



'-n'=0 q=l 



mnq I '^mnq. 




Kanq3n{kqa)P'^kl sin OkdOk 



Bt 



1 



h'h\koa) 



oo 2 



n'=0 q=l 



27iG^nq / A';^^^gj^{kga)P^kl sin ^^d^ 



-Eo[5m,i + Sm,i]b'^Jn{koa) 



(28) 



From those determined scattering coefficients, the radar cross sections (RCSs) of the gy- 
rotropic anisotropic sphere can be calculated, i.e.. 



a = lim 4:7ir 



oo 



•i\2 



4tt 



n=l 



\n ) n 



+ 



US i<l> _ 



sin 9 
B 



A{^e^t> + 



— 7 7* 

n[n + 1) 



nin + 11 



n=l 
Pi " 



1 J n 

de 



At 



sin 9 



^" + 1) 



n[n + 1] 

2" 



(29) 



III. NUMERICAL RESULTS AND DISCUSSION 



To verify this spectral-domain scattering method for the gyrotropic anisotropic sphere, 
we present the bistatic radar cross sections (RCSs) in i?-plane (xoz-plane as shown in Fig. [T]) 
and i7-plane (yoz-plane as shown in Fig. [1]) which are compared to the results calculated 
by a numerical algorithm, i.e., adaptive integral method (AIM)^ extended from Ref[17]. 
The gyro magnetic (e2 = and 7^ in Fig. E^a)) and gyroelectric (62 7^ and yU2 = in 
Fig. ISl^a)) cases have been discussed in Fig. [21 and the good agreement of RCS results on 
both planes is achieved between our method and AIM. It partially verifies that the proposed 
method and the Fortran code developed in this paper are correct. The series in fl26p converge 
rapidly, and it is sufficient to take = 4 as the upper limit of the summation indices n 
and n' . Certainly, it should be pointed out that the convergence rate or the upper limit 
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of the summation depends on the electrical dimension of the sphere (with respect to the 
wavelength) . 




FIG. 2: Radar cross sections (RCSs) versus the scattering angle (in degree) for (a) the gyromagnetic 
sphere and (b) the gyroelectric sphere. The comparisons in RCS results are made between our 
spectral-domain method (solid curve) and the AIM (square dot). The electronic size is fixed at 
koa = vr. 

Then we study a more general case in Fig. |3] in which both material tensors (e and JI) 
are gyrotropic and lossy. The radar cross sections on ii^-plane and iJ-plane have been shown 
in Fig. [3l To the best of our knowledge, the scattering by such a general gyrotropic sphere 
has not been reported, except for its subcase of gyromagnetic spheres^. Obviously, our 
model is more general in terms of the material complexity in Ref[13]. Our spectral-domain 
analysis is distinguished from the spatial-domain method in^^, and one can imagine that 
if the spatial method in Ref[13] is extended to study our general gyrotropic materials, the 
formulation would be lengthy due to the second tensor of permittivity. Hence, even for the 
general gyrotropic materials, our method results in simplified formulation. 

To illustrate the applicability of this analytical solution to the gyrotropic anisotropic 
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FIG. 3: Radar cross sections (RCSs) versus scattering angle (in degrees): The electronic size is 
chosen as A:oa = vr. 





£,=(3+i0.2)8„ 


Mr(2+i0.iK 










B3=(2+i0.1)s„ 


H3=(3+i0.2)n„ 










E Plane ^ 






- - H Plane ^ 
1 . 1 


J k„a=47: 



45 90 135 180 

Scattering Angle (Degree) 



FIG. 4: Radar cross sections (RCSs) versus scattering angle Q (in degree) in i?-plane (solid curve) 
and and -ff-plane (dashed curve). The electric dimension is chosen to be k^a = Att.. 

sphere of electrically large size (for example, in the resonance region), the RCSs of a relatively 
large gyrotropic anisotropic sphere with loss are presented in Fig. 4. The lossy permittivity 
and permeability parameters are chosen as ei = (3 + 0.2i)eo, €2 = eo, es = (2 + 0.1z)eo, 
/ii = (2 + 0.1i)/io, /i2 = /Wo; /^3 = (3 + 0.2i)fiQ. When the dimensions are increased, the 
convergence number (A^ = 24 for the sphere koa = An in Fig. 4) is also increased. 
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IV. CONCLUSIONS 



In this paper, an analytical solution to the scattering by a general gyrotropic anisotropic 
sphere has been obtained. The method is developed based on the multipole expansion of 
the field along with the Fourier transform where the unknown angular spectrum amplitude 
is determined in spectral domain. The three-dimensional electromagnetic scattering of a 
plane wave by an gyrotropic anisotropic sphere has been theoretically formulated, physi- 
cally characterized and numerically discussed. Numerical results for special cases are also 
obtained and verified by comparing with the results from the method of moments. The good 
agreement validates our spectral-domain scattering theory. By using our proposed theory, 
the scattering problems of the general optically anisotropic sphere can be analytically stud- 
ied in spectral domain and RCSs can be readily computed. The analytical solution under 
arbitrary incident angle is still under investigation. 

Appendix A: Scattering coefficients of eigen-expansions in Eqs. (20) and (27) 



Because the spherical wave functions X^„(r, k), iVf^„(r, k), and N^„{r, k) form a com- 
plete set of orthogonal basis functions, we can employ them to expand any solutions uniquely. 



mn 




(A-1) 



e.g., 




Y.[<n%)M^l{r,k,) + hl,^{ek) 



mn 




iV«(r,/c,) + c^J^,)X«(r,/c,)J, 
■ Y.[<n{Ok)M^l{r,k,) + hl,^{e,) 



mn 




N^^^{r,k,) + cl,^{e,)L^^^{r,k,)\, 
■ Y.[<n{Ok)M'^l{r,k,)-rh^^M) 



mn 



N^l{r,K) + ^mn{Ok)L'^l{r,K)\ 



(A-2) 
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The coefficients in (lA-2p . i.e., af^„, and (where p = x,y,z), are functions of 9k 
and (pk- For the detailed expansion and discussion, the information can be found in-^. We 
provide only the coefficients of A^^^, B'^^^ and C^^^ used in the main text. From Eq. (fT5]l . 
we have 



where 



with 



FliOk, = Ff (^,, + FfiOk, (A-3) 



Ff(efc,0,)=F;,^(5,,0fc)x + F,m,0fc)y + F;f(^fc,0,)? (p=l, 2) (A-4) 



-— sm0fc, p = 1, 

— cos0fc, p = 2; 

^COS0fc, p = 1, 

— sm0fc, p = 2; 



qz \ 



k,(Pk) 



0, p=l, 

1, p = 2. 



(A-5) 



In the above equations, the intermediate parameters, Ai, A2 and A, are functions of 
only 6k as given in Eq. ( |T4l) . Then we can split the parameters as follows 

mng mnq ' mnq^ 

Tje rjel i De2 

mnq mnq ' mnqi 

= C^^ -\- C^"^ ( A-(\) 

^mnq ^mnq ~ ^mnqy \ 

and thus obtain (p = 1 or 2) 

Asp -im<f>k _ pep X , pep y . T?^P 
mnq^ qx ^mn ~ qy mn ~ qz ^mn) 

TDep -im<j>k — pephx , jpepuj , pepi^z 
mnq qx mn ~ qy mn ~ qz mm 

riep -im<f>k _ pep x , pep y , pep z ( A-7) 

mnq qx mn ~ qy mn ~ qz mn' \ J 

As a result, we can now obtain the expansion coefficients of E- fields in a gyrotropic 
anisotropic medium, i.e., A'^^^, B!^^^ and C^^g (where q = 1,2), as follows: 
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for p = 1 and m > 



/lel 
mnq 



;el 
mnq 



_ .„ 2n + l (n-m)!Ai 
~ ^ 2n(n+ 1) (n + m)! A 

^-^(cos^,)-P„"^+^(cos^,) 



n-\-m){n — m 



1 



(n — m)! Ai 



"^"5 2n(n + l) (n + m,. _ . 

- ' - --)P™-i(cos^fc) + (n + l 



l)(n + m) 

X (n + m — 1) 
P^V(cos6'fe) + n(n-m+2)(n-m+l) 

X P--/(cos^fe) +nP-V(cos^fe)" 

^ 1 (n-m)! Ai 
2A;g (n + m)! A 



+ m)(n + m — 1) 



X P--/(cosefe) + P„"!.V(cos^fe) 
- (n - m + 2)(n - m + l)P^-/(cos^^fc) 
-P-V(cos^,)' 



(A-8) 



while for p = 1 and m > 0, 



el 



-mnq 



B 



el 

—mnq 



c. 



■el 

—mnq 



(n + m)! 1 
^ (n-m)! 



(-1) 
(-1)"^+ 



m+l + "^J 



(n — mj 
]^ (n + m) 



(n — m) 



TDSl 

mnqi 



c: 



el 

mnq' 



(A-9) 
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Similarly, for p = 2 and m > 0, we have 

.62 ^ n+1 2n + l (n-m)! ;A2 
'"'^^ n(n+l) (n + m)!l2A 



1 m — mi! r 
1) (n + m)! l2A 



+ l)Pr"i(cos^,) + P™+^(cos^fc) 
+ mP„"^(cosefc)}, 



^mnq '' 



n+1 



1 (n — m)! [ Aj 
.2A 



n{n + 1) {n + mj: ^ l 
m-l)P„"r/(cos^ik)-(n + l 



(n + l)(n+m) 



X (n+ 

X P^V(cosejt) +n(n-m+2)(n-m+l) 
X P„T/(cos^,) - nP™V (cos ^,)' 
+ [n(n-m + l)P^i(cos^fe) 

-(n + l)(n + m)P„"r/(cos^fc)]}, 



while for p = 2 and m > 0, 



(^e2 _ -n+i 1 (n-m)! f A^r 

xP™-/(cos^^fc)-P„"!+/(cos^^fc) 
- (n - m + 2)(n - m + l)P;;Xi\cos9k) 
+ P-V(cos^,)" 
xP-(cos^fe)}; 



-1) 



- (2n + l)cos6lfe 



(n-m)!' 

De2 _ / -I \m + pe2 
-mnq I ^ (n-m)! 

^e2 _ / + 



In a procedure similar to the above, the expansion l 
gyrotropic anisotropic medium can be also obtained, 



(A-10) 



(A-11) 



(n + m)! 2 

I mn 

coefficients of the H-field eigenvector in 
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